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Abstract. We attempt to systematically derive perturbative quantum correc- 
tions to the Berry diagonal approximation of the two-level correlation function 
(TLCF) for chaotic systems. To this end, we develop a "weak diagonal approxi- 
mation" based on a recent description of the first weak localisation correction to 
conductance in terms of the Gutzwiller trace formula. This semiclassical method 
is tested by using it to derive the weak localisation corrections to the TLCF for 
a semiclassically disordered system. Unfortunately the method is unable to cor- 
rectly reproduce the "Hikami boxes" (the relatively small regions where classical 
paths are glued together by quantum processes). This results in the method fail- 
ing to reproduce the well known weak localisation expansion. It so happens that 
for the first order correction it merely produces the wrong prefactor. However for 
the second order correction, it is unable to reproduce certain contributions, and 
leads to a result which is of a different form to the standard one. 



1. Introduction. 

Much recent research has concentrated on exploiting the similarities between quantum 
chaos (the quantum behaviour of systems which are classically chaotic) and disordered 
conductors in order to gain a better understanding of both [1-7] For both types of 
system one cannot calculate energy levels and eigenf unctions explicitly since these 
do not have simple analytic forms due to the non-integrability of the systems. 
One therefore concentrates on the statistical properties of the energy levels and 
wavefunctions. The key similarity between quantum chaotic and disordered systems is 
that, in an appropriate regime, their energy level spectra both have statistics described 
by random matrix theory (RMT) [8-10] 

The techniques used to analyse the two systems (quantum chaotic and disordered 
conductors) are very different, and it is this which makes the similarities in the level 
statistics so potentially fruitful. Quantum chaos makes use of the Gutzwiller trace 
formula [11], a semiclassical relation in which all classical periodic orbits are summed 
over to obtain the quantum mechanical density of states. In most work a particular 
system is considered, and statistical averaging is performed over an energy window 
in the spectrum. In disordered conductors, the microscopic details of the disorder 
potential are neither known nor interesting, and this enables one to average over all 
realisation of disorder from the very beginning. This allows one to develop an effective 
field theory which has a well-defined perturbation expansion. As a consequence it is 
possible to calculate statistical quantities in disordered conductors in several important 
regimes. More specifically there are two important length scales in a disordered 
conductor: the system size, L and the elastic transport mean free path, i = vpT, 
where r is the elastic transport scattering time, vp is the Fermi velocity. Using the 
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diffusion constant D — VpT/d, we can generate two time scales: t and the ergodic 
time, Terg — L^/D ~ {L/i)'^T, which is the time taken for an electron to explore the 
whole sample. This leads to the following time regimes: the ergodic regime, r^rg < t, 
where the electron has explored the whole sample; the diffusive regime, Terg < t < t, 
where the electron is diffusing through the sample as it scatters off impurities; and 
the ballistic regime, t < t, where the electron is moving ballistically between scatters. 
These of course lead to related energy regions via the uncertainty relation: ergodic, 
E < Erh] diffusive, Eth < E < fi/r; ballistic, Ti/t < E, where the Thouless energy, 
Exh = ^/^org = fiD/L^. Finally, there is the energy scale A, the mean energy 
level spacing which defines a further energy regime within the ergodic regime: the 
quantum regime, _B < A, corresponding to times longer than the Heisenberg time, 
i > in = ?i/A. In the disordered conductor one can calculate the level statistics 
in all these regimes, [8, 12, 13] and it shows RMT behaviour in the ergodic regime. 
In contrast, the Gutzwiller approach has so far been consistently applied within the 
diagonal approximation developed by Berry [14], although efforts have been made to 
go beyond this approximation by considering action correlations [15]. In the diagonal 
approximation, the leading contribution to the two-level correlation function (TLCF) 
has been found in the ergodic regime of quantum chaotic systems [14], and in the 
ergodic and diffusive regime of disordered systems [1]. Unfortunately, this method only 
works outside the quantum regime; in particular, it does not allow one to calculate 
the corrections in powers of t/tu predicted by RMT in the ergodic regime for systems 
with time-reversal invariance. This might lead one to ask the following questions: (i) 
Are the level statistics in quantum chaotic systems still described by RMT within 
the quantum regime? In particular, do the corrections in powers of t/tu predicted 
by RMT exist in the ergodic regime? (ii) Are there any analogues of the diffusive or 
ballistic regimes in quantum chaotic systems? 

It is the possibility of answering such questions as these that motivates work into 
the analogies between chaotic and disordered systems. 

One approach to answering such questions would be to develop a powerful field 
theoretical technique for quantum chaotic systems analogous to that which exists 
for disordered metals[16, 17]. This field theory should be able to reproduce Berry's 
result as the first term in a perturbation expansion. An effective field theory - the 
"ballistic" sigma-model [2-4] - has recently been developed for quantum chaos. The 
only averaging used was over an energy window in the spectrum, and this seemed 
sufficient to reproduce RMT behaviour. However, this method lacks a well-defined 
perturbation expansion, unless one introduces additional statistical averaging over 
some ensemble [18-20] 

The aim of this paper is to establish whether a regular perturbative expansion can 
be formulated within semiclassical methods based on the Gutzwiller trace formula. In 
particular, we will concentrate on the two-level correlation function (TLCF), R{uj), 
and its Fourier transform, the spectral form factor (SFF), K{t). We consider the model 
of randomly distributed semiclassical scatterers[5, 6]. Standard scaling considerations 
[21] ensure that for sufficiently large systems the results should coincide with those 
derived for a system with Gaussian white- noise disorder [12, 22], which is the 
conventional model for diagrammatic considerations. Note that in the diffusive regime 
in two dimensions the one loop result (corresponding to the diagonal approximation) 
vanishes, and the two loop result is the leading contribution [22] for systems with 
time-reversal invariance, whilst the three loop result is the leading contribution for 
systems without time-reversal invariance. 
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Verifying such an approach for disordered systems should be the first step before 
one can apply the same methods to generic quantum chaotic systems. Perturbation 
theory is used to describe the diffusive regime, whose existence in the disordered 
system is due to the separation of the time scales, r and ierg, by the large parameter 
{L/tf . Such a parameter does not exist in a generic chaotic system. Nevertheless, 
numerically iorg is always much longer than the 'ballistic' time, r. It might be that 
the perturbation expansion of K{t) for t < t < icrg would reveal behaviour which 
is more universal than that characterised just by the contribution from short orbits. 
Perturbation theory is also used to describe the ergodic regime of disordered systems. 
It is known from the exact solution [8] that there are correction in powers of t/tn to 
the leading order diagonal contribution to the SFF (which is linear in t). For such 
systems, perturbation theory is able to reproduce these corrections. However in the 
case of a generic chaotic system, it is not known whether such corrections necessarily 
exist; if perturbation theory could be applied to these systems this could be checked. 

Therefore, we will try to develop a regular perturbation expansion based on the 
Gutzwiller approach. In this case, the SFF is contributed by the sum over all pairs of 
closed orbits which particles traverse in the system. The diagonal approximation [14] 
suggests that only identical (or time-reversed, if allowed) paths make a contribution 
which does not disappear after averaging over some energy interval. The reasoning 
behind this is that different orbits have uncorrelated classical actions, and hence their 
contributions to the SFF have uncorrelated phases and average to zero. 

Perturbation theory for K{t) in the disordered metal [12, 1] gives not only the 
diagonal result in leading order, but also weak localisation corrections [22] at higher 
order. Recently we have shown[23] how to interpret the weak localisation corrections 
for disordered metals in terms of a trajectory picture. It is tempting to use this 
trajectory picture for actual calculations in the language of the trace formula. 

In the trajectory picture, we identified paths contributing to the SFF which are 
identical for most of their length, except at regions where a path self-intersects. At 
these regions it is possible to make two different paths by connecting up the partial 
paths in different ways. We thus obtain a set of contributions to the SFF where 
the actions corresponding to the two paths are almost identical, whilst the paths 
themselves are piecewise identical (see Fig. 1). We introduce the terminology weak 
diagonal approximation to describe such paths in contrast to Berry's strong diagonal 
approximation. The perturbative order of a given contribution is the number of loops 
created by self-intersections (so that the strong diagonal approximation is the one-loop 
term). In disordered metals, the regions where the paths are identical correspond to 
diffusion propagators; the self-intersection regions to the Hikami boxes [24, 17]. The 
latter are regions with size of order £ where s-wave scattering ensures the gluing 
together of diffusive paths. In a general chaotic system, the loops are made by 
gluing together (by diffraction, tunneling or some other quantum process) classical 
trajectories to make a closed orbit. By analogy with the disordered metal, we call 
such a gluing region a Hikami box. We believe that the whole question of whether a 
semiclassical approach can reproduce perturbation theory will reduce to the question 
of whether it can correctly obtain the Hikami boxes - they are the glue which holds 
the coherent partial paths together. 

To carry out the test of semiclassical theory proposed in the last paragraph, we 
need to use a particular semiclassical approach and work up to a particular order of 
perturbation theory. A number of semiclassical methods exist for dealing with what 
we call the Hikami boxes [25-27,5] However since these only calculate the leading 
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order weak localisation correction to conductivity in disordered metals, they deal 
with a relatively simple Hikami box. To test the validity of a semiclassical approach 
one should check whether it reproduces the higher order Hikami boxes. We take 
the approach based on the Gutzwiller trace formula, developed by Argaman[5], as 
such an approach would have a good chance of being directly applicable to chaotic 
systems. We calculate the SFF, and hence the TLCF, to third loop order in the 
perturbation expansion (which involves the calculation of higher order Hikami boxes). 
Note that the two (three) loop order gives the leading contribution to the TLCF in 
the diffusive regime of a disordered system with (without) time-reversal symmetry. 
Unfortunately, this semiclassical method gives the correct functional form only up to 
two loop order. Beyond this, it does not work. In particular, we show that it can not 
reproduce the three-loop result. It appears that this method, although quite general 
and very attractive because of its simplicity and obvious physical interpretation, 
cannot correctly obtain the Hikami boxes. This does not mean that the Hikami boxes 
cannot be reproduced by semiclassical methods: it has been shown [6] that this is 
possible for exactly the system we consider in this paper. However, it is not clear 
whether the technique suggested in [6] may be extended to generic chaotic systems. 
Since the method we consider here is based on the Gutzwiller formula, if it were 
successful, it would be completely general. 

The rest of the paper is organised as follows. In section 2, we discuss the weak 
diagonal approximation and the regimes in which it will apply. In section 3, we 
briefly review the derivation, from the Gutzwiller trace formula, of a semiclassical 
expression for the SFF, and the application of the strong diagonal approximation to 
this expression for the SFF. We then construct the weak diagonal expansion in terms 
of this semiclassical expression for the SFF. In sections 4 and 5 we calculate the two- 
loop and three-loop contributions to the SFF (the first and second terms in the weak 
diagonal expansion) respectively. In sections 6 and 7 we use the results of sections 4 
and 5 to write explicit formulas for the leading contributions to the TLCF of a two- 
dimensional disordered system in the diffusive regime with and without time-reversal 
symmetry. In section 8 we consider whether the techniques used for the diffusive 
regime can be extended to the ergodic regime. Finally in section 9 we discuss our 
results. 

2. The Weak Diagonal Approximation. 

In this section we explain in more detail what is meant by the weak diagonal 
approximation. We then discuss the timescales for which such an approximation 
would be valid. 

The weak diagonal approximation involves evaluating the contribution of paths 
that are nearly-identical. In this context nearly-identical means the pair of paths 
follow each other everywhere, except when the paths come close to themselves. In 
such regions ("Hikami boxes" marked by dashed boxes in Fig. 1) the two paths differ: 
one path crosses in this region, while the other does not. Two paths are considered 
close when their separation in phase space satisfies 5r5p < h. If the system has 
time-reversal symmetry then 6p is either the sum or the difference of the momenta. 

We consider systems in which the Hikami boxes are due to quantum mechanical 
effects such as diffraction, tunnelling, etc. The time scale, tg, of such quantum 
scattering events, will be infinite in the classical limit ?i — > 0, and much larger than 
the flight time (mean free time) r in semiclassical (disordered) systems. Nonetheless, 
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Figure 1. A schematic representation of the first few contributions to the weak 
diagonal approximation of the TLCF. The two paths in each figure follow each 
other exactly, except in the dashed boxes. In the dashed boxes, which correspond 
to Hikami boxes, quantum scattering causes one path to cross while the other 
does not. In (a) the paths are identical everywhere, this gives R\{ll>) which equals 
the strong diagonal approximation. The geometry of paths contributing to R2 (a;) 
are shown in (b). The two different geometries of paths contributing to R3{lo) in 
a system without time reversal symmetry, are shown in (c) and (d). 

The arrows on the retarded Green's functions (R) point in the direction of 
the underlying path, while those on the advanced Green's functions (A) point in 
the opposite direction. 



these quantum events are of great importance for any paths with a period larger than 
tq. These quantum scattering events cause a particle to be spread among a bunch of 
classical paths which are close in phase space. 

For paths with a period much larger than ig , the contribution of a given pair of 
nearly-identical paths will be of the same order as the contribution of identical paths. 
However the probability of a path coming close to itself in phase-space in time t is 
small if i <C in, and hence the weak diagonal approximation provides pcrturbative 
corrections to the strong diagonal approximation. The strong diagonal approximation 
arises from a path which forms a single loop (Fig. la). The leading order correction 
comes from the paths with one self-intersection which form two loops, as shown in Fig. 
lb. The next order correction comes from the paths with two self-intersection which 
form three loops (e.g. Fig. lc,d) and so on. This loop expansion is totally equivalent 
to the standard field-theoretical loop expansion [23]. 

In the ergodic regime, iorg <S i <^ ^h, the probability of a path coming close 
to itself is of order t/tn- Here one expects universal behaviour which is described 
by RMT. The strong diagonal approximation reproduces only the linear behaviour 
of the form factor K{t). The weak diagonal approximation allows the calculation 
of pcrturbative corrections in i/in, thus raising the possibility of checking whether 
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chaotic systems fully obey RMT [8, 9] in this regime. Note that the very possibility 
of this expansion depends on whether the condition tq <C % is satisfied. This is not 
necessarily satisfied for a generic billiard of typical size L, for which tq ~ t{L/Xe)°', 
where the flight time t (x L, and Xe is the De Broglie wavelength of a particle 
with energy E. As in oz L^ for a two-dimensional billiard, the above inequality is 
parametrically valid only for a < 1. 

We have chosen to apply the weak the weak diagonal approximation to a system 
of randomly placed semiclassical scatterers because the above conditions are easily 
satisfied. In this case, the quantum scatterings are mainly due to the paths grazing 
the scatterers, so the theory of penumbra diffraction [28] can be applied. Then for 
a system with N scatterers of radius R per unit volume tq '-^ T{R/XEy^'^, where 
T = {veNR'^~^)~^. This guarantees that tq <C ^h even if the system size L is of 
order R. If we consider an extended system of size L » vet, then tq <^ teig, and the 
weak diagonal approximation can be applied to the diffusive regime of this system. 
In this regime the probability of self- intersection is parametrically smaller than t/tn- 
The quantum parts of trajectories (Hikami boxes) are the same as in the ergodic 
regime whilst the classical parts give a different contribution as they are too short 
to explore all the phase space. When applied to the diffusive regime of a system 
with a Gaussian white-noise disorder potential, the weak diagonal approximation still 
reproduces [23] the standard weak localisation corrections [17, 12]. The system of 
semiclassical scatterers in the limit oi L ^ R must have weak localisation corrections 
of the same form as those for the white-noise disorder; therefore, any consistent 
semiclassical method must reproduce these results. 

It is interesting to ask whether an analogous non-ergodic regime exists for generic 
chaotic systems. For this the time-scales must be arranged as t -^ tq ^ icrg- Since 
T and torg are parametrically the same, the numerical window between them should 
be large enough to squeeze tq inside. It remains to be seen whether such a condition 
may be satisfied by any billiard. 

The final timescale we should consider is the Erhenfest time tsrh — A~^ \n{R/XE), 
where A in the Lyapunov exponent [6]. For a system of semiclassical scatterers it is 
clear that ^Erh ^ tq For a generic billiard, tErh will depend on the logarithm of (L/Xe), 
while it seems likely tq will depend on some power of (L/Xe)- In the semiclassical 
limit, Xe <C R, we expect that such a system will also have r ^ ^Erh <?C tq. This 
will be of relevance to the details of how we construct the Hikami boxes in the next 
section. 

3. A semiclassical description of the strong and weak diagonal 
approximations. 

The two level correlation function (TLCF) is given by, 

i?H = ^Mi?+ic.)z.(i?-i^)\-l (1) 

where t^(-E) is the density of states per unit volume. Given the spectrum of the system, 
{E„}, v{E) can be written as, 

^{E) = ^^Y.^{E-E^) (2) 
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and the average density of states, v = {v{E)) = (AL'') . The averaging denoted by 
(• • •) is carried out over a certain energy window in the spectrum of the chaotic system's 
spectrum. In disordered systems it is also possible to average over the ensemble of all 
realisations of the disorder. Ensemble averaging is used in the vast majority of work 
on disordered systems, because it is simpler and better understood. In both chaotic 
and disordered systems, the mean level spacing is given by the Weyl rule, so a system 
of volume, L"*, at an energy, E, has A = h'^/riE- The volume of the phase space 
constant energy surface is Qe — {SdPE^'^L'^)/vE, where Sd is the surface area of a 
d-dimensional sphere of unit radius, and pE and ve are the momentum and velocity 
of a particle with energy, E. 

We consider a system of non-interacting electrons in a potential which varies 
slowly on the scale of the wavelength of the electrons at the Fermi surface. In this 
semiclassical limit, the density of states is given by the Gutzwiller trace formula [11], 

v{E) -v = Y, (mE) exp [iS^{E)/h] - Al{E) exp [~iS^{E)/h] \ (3) 

where the summation is over all periodic classical paths, Aa{E) and Sa{E) arc the 
amplitude and action of the ath classical path. For convenience we have defined 
Aa{E) so that it includes the phase factor due to the Maslov index. When studying 
spectral statistics semiclassically, it is often more convenient to consider the spectral 
form factor (SFF), K{t), than the TLCF itself, 

K{t) - ( ^ A^Al exp [z(5„ - Sf,)lh] 5[t-\{T„+ T^)] \ (4) 

where the period of the classical periodic path, a, is Ta — dSa{E)/dE. The TLCF is 
related to the SFF in the hmit that uj <^ E,hy 



R{uj) = 4A2sRe 



dt exp [iojt/h] K{t) 



(5) 



Berry introduced the diagonal approximation of (4) [14], arguing that after the 
averaging was carried out, only terms where a — (3 would still contribute to the 
summation. All terms with Sa ^ Sp oscillate wildly as the averaging is carried out, 
and so these terms will be negligible. Berry's diagonal approximation assumes that 
those terms for which a ^ [3 but Sa ~ Sp are rare enough that one can ignore them. 
We will call this approximation the strong diagonal approximation, to distinguish it 
from the weak diagonal approximation that we will introduce below. When applied to 
the ergodic regime in a generic chaotic system, the strong diagonal approximation was 
found to reproduce the leading order behaviour of the spectral form factor in random 
matrix theory. More generally, it was shown that the contribution of these strong 
diagonal terms to the SFF is given by [1], 

where /? = 1 if the system has time reversal symmetry, and /? = 2 if it does not. 
Here P{t) is the probability of returning to the same point in phase space, in a time, 
t, integrated over the phase space surface of constant energy, E. We now introduce 
notation that will be of use later on, let /^(r',^', t; r,p) be the probability that a 



Can the trace formula describe weak localisation? 8 

particle which is initiaUy at (r,p), is at (r',p') after a time, i, given that the particle's 
final energy, E^ is constrained to equal its initial energy. Then (6) becomes 

2 t f 
Ki{t) = - /2^^p / drdpE fs {r,p,t;r,p) (7) 

where dps = dpd [H{r,p) — E], in other words the integral is restricted to the surface 
in phase space with energy, E. When applied to the diffusive regime of a two- 
dimensional disordered system, this result is t-independent, which means it gives no 
contribution to the TLCF. 

We now attempt to develop the weak diagonal approximation for a semiclassical 
system based on the Gutzwiller Trace formula. To do this we must consider how to 
deal with the regions where the paths are not identical (Hikami boxes) . A path a that 
comes close to itself, will typically stay close to itself for a time of the order of the 
Erhenfest time ^Erh- Since ^Erh ^ ^q, the nearest quantum scatterer on the path will 
be some distance from the region where the path comes close to itself. Therefore the 
behaviour of path a in the Hikami box can be considered as classical. Path f3 is the 
path that will cross in the Hikami box, it will be identical to a up to the last quantum 
scattering before the region in which a is close to itself. Path /3 must leave that 
scatterer with a slightly different momentum from path a, and so moves away from 
it slowly before converging slowly towards the other leg of path a. Finally /3 becomes 
close enough to a that at the next quantum scattering the two paths become identical 
again. When path a comes close to itself, the distance between the two parts of the 
path is of order of a wavelength, A. Path /3 is never more than that distance away 
from one or the other of the legs of path a. These distances are much smaller than the 
lengthscale of the potential, R. Therefore the scatterings of path /3 are all correlated 
to scatterings on one or other leg of path a. This means that ensemble averaging 
will not destroy the contribution of this pair of paths. To evaluate the contribution 
of pairs of paths of this type we construct a method based on that used by Argaman 
[5]. Since we are assuming that for every path a, there is a path /3, where a and (3 
are described above, we can assume the amplitude of path (3 is approximately equal 
to that of path a. However path /3 has to go slightly further than path a, because it 
crosses itself. This means we expect the action of path (3 to be slightly greater than 
that of path a [5] . 

Note that the Hikami boxes in the standard diagrammatic approach are more 
complicated: they are 'dressed' by additional impurity lines, which is equivalent to 
having an extra quantum scattering inside the box. It is not clear why the contribution 
of such an extra scatterer should be important in the semiclassical approach, since its 
inclusion would be a small correction to the semiclassical Hikami box described above. 
However, the two methods are so different that such a direct comparison does not make 
much sense. 

We now make the weak diagonal approximation, by assuming we can ignore all 
terms in the double summation in (4) which do not have the geometry described above. 
Clearly terms in the double sum with f3 ~ a have already been taken into account by 
the strong diagonal approximation. We can then substitute the following into (4), 

Ap ^A^ 5/3 = 5„ + SS (8) 

Since we are considering pairs of paths that are not identical, the double summation 
over a and (3 in (4) gives two contributions for each pair of paths. For example if 
we label a given pair of paths, 1 and 2, then one gets a contribution to the double 
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Figure 2. The Hikami boxes (the dashed boxes in Fig. 1) are shown here in more 
detail. In (a) ri and r2 are close, while pi and p2 are approximately anti-parallel. 
In (b) r\, V2 and r^ are close, while pi, p2 and pa are approximately parallel. 
The directions of the arrows are explained in the caption to Fig. 1. 



summation for a = 1 and (3 — 2, and an identical contribution from a = 2 and P — 1. 
The assumptions in (8) allow us to write the ith term in the weak diagonal expansion 
of the spectral form factor as, 



K,{t)^2y(\A^fS{t-Ta,Jcxp[~iSS^/h] 



(9) 



where the sum over a, is over all periodic paths have i loops. As can be seen in Fig. 1, 
to derive K2 (t) and K3 (t) we need the action difference in two geometries of Hikami 
box, these geometries are shown in Fig. 2. The first is one in which each path comes 
close to itself once (Fig. 2a) and the second is one in which each path comes close to 
a single point on itself twice (Fig. 2b). In the first case, each path enters the Hikami 
box twice, and the action difference between the two paths is, 

SS2 = {ri-r2)-{p2+Pi) (10) 

In the second case, each path enters the Hikami box three times, the action difference 
is, 



SS3a = {ri - r3).pi + (r2 - ri).p2 + (^"3 - r2).p3 



(11) 



The reason for the notation we have chosen for the subscripts on 6S2 and SS^a will be 
made clear later in this paper. 



4. The two loop contribution in the weak diagonal expansion. 

Here we calculate two loop term in the weak diagonal expansion. This is the leading 
order correction to the strong diagonal approximation (1 loop) result. The geometry 
of paths that contribute to the two loop term is shown schematically in Fig. lb. 
We expect this to be the dominant contribution to the weak diagonal perturbation 
expansion in a system with time-reversal symmetry. However in a system without 
time-reversal symmetry the two loop term in the weak diagonal expansion is zero. 
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This is because one path follows the time-reverse of the other path for some of its 
period. In a system without time-reversal symmetry it is impossible for a path to 
exactly follow the time-reverse of another path. 

In all the figures in this paper we have followed the Feynman diagram convention 
for the arrows. This means that the arrows on the retarded Green's functions (R) point 
in the direction of the underlying classical path, while the arrows on the advanced 
Green's functions (A) point in the opposite direction. This leads to the slightly 
counter intuitive observation that those pairs of paths which have the arrows pointing 
in opposite directions at all points are the only ones to survive in a system without 
time-reversal symmetry. In the language of disordered systems, the pairs of paths 
which have the arrows pointing in different directions are diffusons, while those with 
the arrows pointing in the same direction are Cooperons. 

The two loop term in the weak diagonal expansion, shown schematically in Fig. 
lb. can be found by substituting (10) into (9), 



K2{t) - 2^( |A„|' 5 {t - Tc.) cxp 



1 

--(ri -ra) • (P2+P1) 



(12) 



where the ath path is a path which starts at (r,p), passes through (ri,pi) at time 
ii, and (7*2,^2) at time ^2, before returning to (r,p) at time, t. The sum over a is 
over all primitive periodic paths and all repetitions of them, however in what follows 
we will ignore the possibility of repetitions. This is justified in [1] by noting that the 
number of primitive periodic paths grows exponentially with the period of the path. 
Therefore repetitions make a very small contribution to the sum over all paths with a 
given period. 

The amplitude A^ is the same as in the derivation of Ki{E, t), the following result 
was derived in [1], 

\A^\'' 5{t-T^) ^ —^ J drdpE fMr,P,t;r,p) (13) 

where /^^(r', p', t; r, p) is the probability of a particle following path a from (r, p) to 
{r',p') in a time t, given that the particle's final energy, E, is constrained to equal its 
initial energy. The integral over momentum is restricted to the energy surface, dp^; is 
defined below (7). Equation (13) was originally derived for purely classical paths. We 
want to consider paths which arc mainly classical, but which occasionally undergo a 
quantum scattering. For purely classical paths /^^(r', p', t; r, p) is simply a (5-function 
at the point in phase space which a classical path starting at (r,p) would reach after 
time, t. We include the occasional quantum scattering event by dividing a given path 
up into segments of classical and quantum behaviour. The classical segments will 
provide most of the path, but there will be occasional short quantum segments. The 
classical segments can be related to the classical probability of propagating between 
two points in a given time. The quantum segments of the path have to be dealt with 
by finding the quantum mechanical propagator for that scattering event. If a quantum 
segment of the path goes from (ri,pi) to (ri+i,pi+i), this can be written in the form 
j4exp[«i?] where A and R are defined by the quantum mechanics of the scattering 
event, and are functions of (rj^i+i,pj^i+i). An example of such a quantum scattering 
event is the penumbra diffraction of paths that graze semiclassically large scatterers, 
discussed in [28]. The probability of propagating from {vi^pi) to (ri+i,pi_|_i) is simply 
\A\'^ . Since the scattering is not classical, this probability will not be a (5-function. 
So a particle on a classical path encountering this quantum scatterer will be spread 
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over a number of classical paths after leaving the quantum scattering region. The 
probabilities of propagating along each segment of the path can be multiplied together 
to give the probability of following the whole of that path. This procedure gives a 
result of the same form as (13), however the probability /^^(r'jp', t; r,p) can now be 
written down for any path, whether it is purely classical or includes quantum scattering 
events, such as diffraction, tunnelling, etc. Similarly, when the sum over a is carried 
out it is now a sum over all the paths including quantum scattering events rather than 
just purely classical paths. In the limit of paths much longer than ig, a particle can 
be considered as mostly following classical paths, but with a small amount of diffusion 
from one classical path to those nearby in phase-space [6] , where the timescale for this 
diffusion is set by ig. 

Consider a quantity which is only non-zero in the vicinity of the ath isolated 
periodic path through the system. If one ignores repetitions of the path then the 
integration of that quantity over the surface in phase space with energy, E, can be 
written in terms of coordinates parallel, (r||,p||), and perpendicular, (rx,px), to the 
path. Then, 

/"dr (•••) = fdr^^dr^ {■■■)= VE f dt' fdr^i---) 

/ dpE (•••)= / dp||dp_L Vj^^6 [p|| - pe] (•••)= Vj^^ / dp± (■ • •) 

where the dpE is defined just after equation (7). This was used in [1] to prove (6), 
and will be of great use to us. 

The exponent in (12) is in terms of the two points (t"i,Pi) and (r2,P2), which are 
both somewhere on path a. It is therefore necessary cast the right-hand side of (13) 
in a form which specifies that these points are indeed on path a. We will now discuss 
the details of how this is done. Consider a path which starts at (r,p). We want to 
know the probability that the path passes through {r',p'). We do not care when the 
path passes through {r',p'), so long as it does so at some time between time ta and 
time tb- So long as the period of the path is more than t2, the probability density 
for this to occur is J^'' dt' faE (f' , p' , t' ': f , P) ■ Therefore the probability of going from 
(ro,po) to (r,p) in time, t, and passing through {r',p') on the way from {ro,po) to 
(r,p) is, 



(14) 



/ di' JaE {r,p,t - t' ]r' ,p') JaE {r' ,p' ,t' ]rf),pf)) 
Jo 



(15) 



Given that the propagation probability for the path a from (t"o,Po) to {r,p) in time i, 
is faE {i^,P,t]ro,Po), one can insist that the point {r',p') is somewhere on the path. 
Then casting all coordinates in terms of components parallel and perpendicular to the 
path, 

faE{r,p,t;ro,po) = / dr^dp'j^ 

X / dt'faE{r,p,t~t';r\p')faE{r\p\t';ro,Po) (16) 

The integral over all (r^,pY) in the vicinity of the path a will pick up a 5- function 
when {r' ,p') is on the path, with no contribution the rest of the time. 
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As previously mentioned we wish to specify that the points {ri^pi) and (^"2,^2) 
are somewhere on path a. Using (16) to specify this, (13) becomes, 

\Aa\^ 5{t- Ta) = / drdpE / dri^dpi^ dr2±dp2± 

X dti dt2 .faEir,p,t;r2,P2,t2;ri,pi,ti;r,p) (17) 

Jo Jti 

where the integrals over ri.2 and Pn^2)E ^i'g over the phase space surface with energy, 

E. We have defined faE {'r',p,t;r2,P2,t2]ri,pi,ti;r,p) as the probability that a 

classical particle with energy, E, which starts at (r,p), passes through (ri,pi) at 

time, ii, then through {r2,P2) at time, i2j before returning to {r,p) at time, t. Hence 

faE{r,p,t;r2,P2,t2;ri,pi,ti;r,p) 

= faE {r,p,t - t2;r2,P2) faE {r2,P2,t2 - ti;ri,pi) JaE {ri,pi,ti;r ,p) (18) 

The integrals over r and pe are carried out by noting that, 

/ drdpE faE {r,p,t;r2,P2,t2;ri,pi,ti;r,p) = f^E {ri,pi,t;r2,P2,t';ri,pi) (19) 

where we have defined t' = t2 — ti. The integrand is now independent of ti, so the 
integral over it simply generates a factor of i. Substituting (19) and (17) into (12), we 
note that we only expect a contribution when p2 is approximately anti-parallel to pi . 
Therefore we can resolve r2 and p2 in the exponent of (12) into components parallel 
and perpendicular to the momentum pi. The integrals over r2j_ and p2i_ can then be 
evaluated using the following stationary phase approximation, 

dxdpf{x,p)exp[ixp/h] ^ (27rn)/(0,0) (20) 

This means the the only contribution comes from paths where r2j_ ~ and p2i_ = 0. 
The integral over P2|| is constrained by the energy (5-function, so \pi\ = \p2\ = pe, 
and therefore P2 — Pi ■ Hence we find, 

\Aa\'^S{t^Ta)ex:p[-iSS2/h] 

= (2^fi)3-d / dT-i_Ldpi_L / dt'faE{ri,Pi,t;ri,-pi,t';ri,pi) (21) 

where SS2 is given by (10). The integral over rix and pix can be converted back to 
the system's coordinates by noting that since the integral over t' makes the right-hand 
side of (21) independent of rf, , the result is unchanged if we integrate over r',, , so long 
as we divide by the path length, VEt. Then (14) can be used to put the integrals in 
terms of the systems' coordinates. However when we do this, the integral of ri and 
PiE over all phase space double-counts all the paths, so we have to divide through by 
two. Since we get a double contribution from each path if the system has time-reversal 
symmetry and no contribution if it docs not, we multiply through by 2(2 — (3). Hence, 

\A^f Sit- T^)cxp [-1682/ h] 



(2 - P)t <■ '•* 

(27r?i)3-- J J, 



dridpiE / dt'faE{ri,Pi,t;ri,~pi,t';ri,pi) (22) 

Finally substituting this into (12), summing over all possible paths, and carrying out 
the averaging (• • •), we find, 

K2{E,t) ^'^^—^ J dridpiE J dt'(/£(ri,pi,i;ri,-pi,i';ri,pi)) (23) 
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At present we have not specified what form the averaging (• • •) will take in the 
definition of the SFF and the TLCF. However it is clear from (23) that it only affects 
the propagation probability. So if the averaging were over an energy window near the 
energy, E^ then one only needs to know what the propagation probability is after it 
has been averaged it over that energy window. In this paper we will be applying the 
above result to a disordered system, so we can chose to define the averaging in the 
definition of the TLCF (1), and hence the SFF (4), as averaging over an ensemble of 
systems with different positions of the scatterers. In this case the averaging of the 
propagation probability in (23) would be over this ensemble of systems. 

Equation (23) is the two loop term in the weak diagonal expansion of the SFF. 
When we set /3 = 1, it applies to a system with time-reversal symmetry, and is the 
leading order correction to the strong diagonal approximation. When we set /3 = 2, 
it applies to a system without time-reversal symmetry, for which the two loop term is 
zero. Therefore in the (3 = 2 case, it is the three loop term which is the leading order 
correction to the strong diagonal approximation. 

5. The three loop contribution to the weak diagonal expansion (for a 
system without time- reversal symmetry). 

In this paper we are interested in calculating the leading order correction to the 
strong diagonal approximation. As discussed at the beginning of the previous section, 
the two loop term of the weak diagonal perturbation expansion is zero in a system 
without time-reversal symmetry. Therefore, to find the leading correction to the strong 
diagonal approximation for a system without time-reversal symmetry, it is necessary 
to calculate the three loop term in the expansion. In general there are five path 
geometries that contribute to the three loop term in the expansion, however we are 
only interested in the three loop term for a system without time-reversal symmetry. 
In such a system three of these five possible geometries give no contribution because 
one path follows the time reverse of the other at some point in their period. This 
means we have the contribution of two geometries to calculate, these geometries are 
shown in Fig. lc,d. 

The first contribution is shown in Fig. Ic, the action difference between the two 
paths, (55*30, is given by (11). The amphtude, A^, is given by taking (13) and using 
(16) to insist that (ri,pi), (r2,p2) and {r^,p^) are somewhere on the path a, 

\Aa_\^ 5{t- To) = I ArdpE / dri^dpix Ar2^dp2i_ dra^Ldps^ 

X / dti / di2 / dts faEir,p,t;r3,p3,t3;r2,P2,t2;ri,pi,ti;r,p) (24) 

Jo Jti Jt2 

The integral over (r,p) is carried out using (19). We only expect contributions when 
P2 and P3 are approximately parallel to pi. So substituting (11) and (24) into (9), 
we carry out the stationary phase integrals over (r2_L,P2_L) and {r3±,p3±), this leaves 
the integral over (ri_L,pi_L). As in the two loop case we introduce an extra integral 
over ri\\ and divide through by the length of the path, w^t, the integral of ri and piE 
is then over the phase space surface with energy, E. However now this integral over 
all phase space counts each contributing path three times, so we must divide through 
by three. For completeness we also multiply the result by the time-reversal symmetry 
factor, 2//3, so that when /3 = 2 (/3 = 1) we count each path once (twice). However 
one should note that to calculate the total contribution in a system with time-reversal 
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synimetry, one would not only have to set /9 = 1 in these equations. One also has to 
calculate the contribution of the trajectories that we ignore because they are zero in 
systems without time-reversal symmetry. Summing over all paths and carrying out 
the averaging (• • •), we find, 

At f /■* f*"*! 

^-(^' ^) - m2^W^^ J ^^^^^- 1 < I '''^ 

jE(ri,Pi,t;ri,pi,t2;ri,pi,t[;ri,pi)) (25) 



where we have defined t'l = t2 — ti and t'2 = ^3 — ^2- 

The second contribution to the three loop term is shown in Fig. Id, the action 
difference between the two paths is, 

(iSsb = (ri - rg) • (p3 - pi) - (r2 - ri) ■ (p4 - pa) (26) 

and the amplitude can be written in terms of probabilities, by taking (13) and using 
(16) to insist that {ri,pi), (r2,P2), ('"3,^3) and (r4,p4) are somewhere on the path 
a, 

|Aq|^ 6{t- To) = . .^ / drdpE / dri^dpi±_ dr2±dp2i_ dr3^dp3±_ dr4^dp4^ 

rt rt+ti j-t+ti I't+ti 

X / dti / di2 / di3 / di4 

JQ Jti Jt2 Jts 

>^.faE {r,p, t; r4, p4, tr, r^,p^, t^; ra, P2, ^2; ri,pi,ti;r,p) (27) 

The integral over (r, p) is carried out using (19). We only expect contributions when p3 
is approximately parallel to pi , while p4 is approximately parallel to p2 ■ Substituting 
(26) and (27) into (9), we carry out the stationary phase integrals over (^34^,^3^) and 
(r4,p4). This leaves the integral over {ri±,pi±) and (r2_L,P2±), wc turn these into 
integrals over the phase space surface with energy, E, by introducing extra integrals 
over riH and r2\i, while dividing through twice by the path length. However now the 
integrals of ri.2 and Pn^2)E over all phase space count each contributing path four 
times, so we must divide through by four. As before, we also multiply the result by 
the time-reversal symmetry factor, 2//3. Summing over all paths and carrying out the 
averaging (•••), 

K^biE, t) - p,^\i.2d J dridpiB dr2dp2B J dt[ J ' dt'2 j ' ' d^ 

/ij(ri,pi,i;r2,P2,t3;n>Pi>*2;^2,P2,ii;ri,pi)) (28) 



where we have defined t'l = t2 — ti, t'2 = t^ — t2 and t'^ = t4 — t^. The sum of (25) 
and (28) gives the three loop term in the weak diagonal perturbation expansion of the 
SFF for a system without time reversal symmetry. Since for such a system the two 
loop term in the expansion is zero, this is the leading order correction to the strong 
diagonal approximation. 

6. The leading order behaviour of the TLCF in a two-dimensional 
disordered system with time-reversal symmetry. 

We will now apply the results of the weak diagonal expansion to a two-dimensional 
dilute system of randomly placed semiclassical scatterers. Each scatterer is 
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semiclassical in the sense that its radius, a, is much greater than the wavelength, 
\p. The typical distance between the scatterers is much larger than their radius. 

In any given system the scatterers are randomly placed. We will be ensemble 
averaging our results, where the ensemble to be averaged over is an ensemble of 
systems with the same macroscopic properties but with different positions of the 
scatterers. The propagation probabilities through two separate regions of the potential 
are uncorrelated because the scatterers are randomly placed, and can therefore be 
ensemble averaged separately. We will assume, 

{Je iri,Pi,t; ri, -pi, t'; ri,pi)) 

-(,/B(n,Pi,t-i';r-i,-Pi)) {fE{ri,-pi,t';ri,pi)) (29) 

where (• • •) now denotes ensemble averaging. This assumption will be valid so long as 
t,t' ^ isrh, because then the region over which the parts of the paths are correlated 
is a small proportion of the whole path. 

On these relatively long timescales, the classical behaviour of the system is 
diffusive. The propagation probability, /£;(r',p',i; r,p), is not entirely classical, it 
include some diffractive scatterings, however these are not going to change the diffusive 
nature of the propagation probability distribution. After a small number of scatterings 
the direction of the momentum of a classical path is effectively randomised. Therefore 
the ensemble averaged propagation probability is independent of both the initial and 
final direction of the momentum. Hence, 

{fE{r',p',t;r,p)) =W{r\r;t)-^ (30) 

^dPE 

where S^ is the surface area of a d-dimensional unit sphere. W{r,r';t) is the 
probability of a diffusing particle which starts at r being within dr' of r' after time, 
t, it therefore satisfies the diffusion equation, 

§i-DWl, W(r\r;t)^6{t)6{r'-r) (31) 

For time scales smaller than the ergodic time, torg — L^/D, the system is in the 
diffusive regime, the diffusing probability distribution W{r', r; t) does not know about 
the boundaries of the container and, 

Wir',r;t) ^W{r'-r;t) = {iTrDty^/^ cxp \-\r'-rf/{4:Dt)\ (32) 

Disordered systems are believed to obey one-parameter scaling, hence it is usual to 
characterise disordered systems in terms of a single parameter, their dimensionless 
conductance, g. This is simply the conductance measured in units of e"^ /(udh), where 
the parameter Ud — tt{2ttY /Sd- The dimensionless conductance of a disordered system 
with volume, L'^, and diffusion coefficient, D, is given by, 

g = UdEThl^ (33) 

where the Thouless energy, Exh = S-Z^org- We will therefore write our results for 
disordered systems in terms of the dimensionless conductance. 

If one uses this description of a diffusive system to calculate the strong diagonal 
approximation of the SFF (6), one find that in two dimensions the result is t 
independent. Since the TLCF is the Fourier transform of the SFF, the strong diagonal 
approximation of the TLCF is zero. This means that it is the corrections to the strong 
diagonal approximation that give the behaviour of the TLCF in a two dimensional 
diffusive systems. Here we will apply the weak diagonal approximation results derived 
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in the previous sections to a diffusive system with time-reversal symmetry. By using 
the above description of propagation in a diffusive system, we calculate the two loop 
term in the weak diagonal perturbation expansion for such a system. This will give 
the leading order behaviour of the TLCF in a two dimensional diffusive system with 
time-reversal symmetry. 

We substitute (30) and (23) into (5), and find the semiclassical evaluation of the 
two loop contribution to the TLCF, R^[uj), is given by, 



Rfiio) = 8(2 - P)L' 



»e 



\2TThJ 



/ dtt exp[iujt/n] / dt'W{0-t')W{0:t~t') 
Jo Jo 



(34) 



where the level spacing, A, is simply given by A = h'^'VE/ {SdP% ^L"^), where 
ScLPe'^L'^ /ve is the volume of the constant energy surface of the phase space of the 
system. 

Thus substituting (32) into (34), gives Rfi^) ^^ the diffusive regime, 

8(2 -/3)£^'^ /Ay 



i?r(u;) 



Ke 



/ Att exp[iwt/n] / At't'-^/'^{t-t'y'^''^ 
Jo Jo 



(35) 



One can evaluate the integral over t' by noting it is a Eulcr /3-function and so can be 
written in terms of F-functions. The result in 2 dimensions can be found by setting 
d = 2 + e, and expanding the F-functions in e. 

However to compare the results of the semiclassics with those of the impurity 
diagram technique, we can Fourier transform the 'Ws in (34) and then compare 
directly the functional forms of the results in the two approaches. If wc define W {q, oj) 
using, 

^'^', . ,,,, 

,c^) (36) 



/AdQ Af. 
, ;.+, cxp [-1 {q-r + Lot) /h] W {q, 



{2TTh) 

then in the diffusive regime, 

W{q,Lu)^h^[Dq^ -ihuj]'^ (37) 

Substituting (36) and (37) into (34), and carrying out the integrals over t and t' , we 
find the semiclassical result for the two loop term in the weak diagonal perturbation 
expansion. 



Rl^{u;) = 8{2-f3)hL'''f^) 3?e 



d{iuj) 



(38) 



in deriving this we have defined, 

d'^q 



Ini^o) 



n^" \Dq^ 



ihuj] 



{2T:hy 
F(n-rf/2) fnYfn.y/^ 



-i) 



d/2-r. 



(39) 



(47r)d/2(n-l)! \aJ \L^gJ 

where g is the dimensionless conductance given by (33), and A is the average level 
spacing. 
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Figure 3. This dictionary allowing a translation from the trajectory pictures 
to the diagrams in the impurity diagram technique. A region where the two 
paths are identical corresponds to a diffuson (or a Cooperon, if one of the arrows 
is reversed). The semiclassical Hikami boxes arc shown with the corresponding 
Hikami boxes in the impurity diagram technique. 



We will now calculate the equivalent term in the impurity diagram technique 
[29, 17, 24]. The trajectory pictures used in the semiclassical approach (see Fig. 1 
and 2) can easily be mapped onto the diagrams in the impurity diagram technique 
[23]. A dictionary for converting between the two is shown in Fig. 3. However this 
correspondence is fairly superficial, the picture drawn in the two approaches maybe 
similar, but the approach to calculations is different. Two dimensions is the critical 
dimensionality for the impurity diagram technique [16]. As a result rcnormalisation 
ideas apply, and one can ignore short trajectories and short range behaviour (in other 
words the large uj and q behaviour). By ignoring these ultraviolet divergences, the 
impurity diagram technique calculations are greatly simplified. 

The two loop term in the impurity diagram technique is given by the trajectories 
shown in Fig. lb, which can be converted to the more usual diagrammatics picture 
using Fig. 3. After throwing away the ultraviolet divergences, the two loop 
contribution is, 



i?2(w) =A{2- l3)nL^ 



\2TTn 

carrying out one of the differentials over lu, 



5fte 



d^ 



R2{uj) =4(2-/3)?ir 



f— V 

\2'KhJ 



sRe 



d 



d{iuj) 



_d{iujf 



2iuj _ 



\lL0{h{uj)f 



{h{uj)r + ^hHhiio) 



(40) 



(41) 



Since the divergent terms can only be ignored in the vicinity of the critical dimension, 
d = 2, one must remember that (40) and hence (41), are only true for d — 2 + e where 
e is small. 
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The functional form of the resuhs of the two approaches, (38) and (41), are 
different. The impurity diagram technique result has an extra term, however it 
so happens that in two dimensions, this term is purely imaginary and so does not 
contribute to the TLCF. This is most easily shown by evaluating its contribution in 
d ~ 2 + e, writing the integrals in terms of F-functions and then taking e to zero. 
Therefore in two dimensions the two approaches have the same functional form. In 
two dimensions the semiclassical result simply disagrees with the impurity diagram 
result by a factor of two. Since the impurity diagram technique result (41) is only 
true for two dimensions, our conclusion is that the semiclassics gets the right answer 
but with the prcfactor wrong by a factor of two. However it looks suspiciously like it 
is so close to the correct result only as the result of a lucky cancellation of the term in 
the impurity diagram technique that the semiclassics can not reproduce. To confirm 
if this is the case we need to go to the next order in the weak diagonal expansion, and 
compare that result to the three loop term in the impurity diagram technique. 

For the sake of completeness we now evaluate R^{lo) in d = 2, taking (38) setting 
d — 2 + e where e is small, 

(i^l«. [,_„,.-■] ^ -<|^ 

7r(^^+e £ L J 2g \s\ 

where Iiiuo) was evaluated using (39), g is the dimcnsionlcss conductance given by 
(33), s ~ w/A, /? = 1 (/3 = 2) for a system with (without) time- reversal symmetry. 
As we have stated, this result is a factor of two larger than the equivalent impurity 
diagram result. 



Rf{s^) = - 



(42) 



7. The leading order behaviour of the TLCF in a two-dimensional 
disordered system without time-reversal symmetry. 



We now evaluate the first of the two semiclassical contributions to the three loop term 
in the weak diagonal expansion of the TLCF, i?3^(w), for a disordered system without 
time- reversal symmetry in the diffusive regime. Substituting (29) and (30) into (25), 
we find. 



RU^) 



3/3 V2^ 



Ke 



dt t exp [iujt/1i\ 



X / dt'idi2W^(0;i'i)T4^(0;4)T4^(0;i-t'i-t2) 
/o 



(43) 



In the diffusive regime, one can use (32) to show this becomes, 

-3d / A \ 4 

Rfaioj) 



16L3 



(. 



A 



3/3(47ri:»)3^/2 \^27rn 



5Re 



dt t exp [iujt/tt\ 







f* r 1 -d/2 

X / dt'^dt'^ t\t'2 {t-t'2- t'l) 
'- 



(44) 



However as before we Fourier transform (43) to the (q, cij)-representation to allow 
comparison with the impurity diagram technique. Substituting (36) and (37) into 
(43) and carrying out the integrals over t and t'-^ 21 



RVai^) 



3/3 



A 

271% 



Ke 



d 

d{iijj' 



[h{u^)r 



(45) 
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In (2 + e) dimensions where e is small, this becomes, 



i?3a(sA) 



(46) 



We now evaluate the second of the two semiclassical contributions, R^Klj), 
for a disordered system without time-reversal symmetry in the diffusive regime. 
Substituting (29) and (30) into (28), we find. 



Rllio^) 



AL^"^ ( A 



27r?i 



/3 



5Re 



dr 



di t cxp \iLotlf\\ 



At\ 



t-t' 



dtn 






dt'o 



xW{r;t[)W{-r;t'^)W{r;t':,)W{-r;t-t[-t'^-t'3) 
where r = r2 — ri. Substituting in (32) this becomes. 



(47) 



Rlli^) 



/3(47rL>)2d 1^2^ j 



di'i 



5Rc 



dr 



di i exp [iciji/?i] 



^ dt' 











dt'. 







X exp 



4L> vt'i t; 



^1^2^3 (^ ^3 ^2 ^l) 
1 



-(i/2 



- ' _,_,_, , (48) 

However we proceed by Fourier transforming (47) to the (q, w)-representation to allow 
direct comparison with the impurity diagram technique. Substituting (36) into (47) 
and carrying out the integrals over r, t and t'l 2 3- 



Rlliuj) ex sRe 



,i7r/2 



3d/2-4 



/"d'^q-id'^q2d'^q3 



{ql + 1) (q2 + 1) (ql + l) f(q, + q^ + q,f + 1 



(49) 



For d = 2, the integral is convergent, so the term in the square bracket has no real 
part, and hence R^'j^{uj) ~ 0. 

Since the semiclassical result for /Jg^ is zero, the total semiclassical result for the 
three loop term of the expansion, in a (2 + e)-dimensional disordered system without 
time-reversal symmetry, is given by (46) with (3 — 2. Expanding the in powers of e, 
and taking the limit e ^ 0, we neglect the unphysical 1/e divergent term. Then in 
2-dimensions, 

«f <'^' = w-w <=»' 

where s = uj/A. We now compare this result with the impurity diagram result for 
i?3(a>), in a system without time reversal symmetry {(3 = 2). The geometries of paths 
that contribute to i?3(cij) are shown in Fig. lc,d, these can be put in terms of the usual 
impurity diagrams using Fig. 3. In dimension d = 2 + e where e is small, one finds. 



i?3(sA) = -TiL 



3d 



f-V 



ih 



b') 



5Re 



92 



_d{i(jj)'' 



[lL0{h{iu)f 



(51) 
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Using the result for /i(cj) given in (39), this becomes, 
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This resuh has a (— |s|~^) dependence, while the semiclassical result has a (In |s|/|s|) 
dependence. Therefore the semiclassical approach gives a result which is at odds with 
the impurity diagram result. 

We find it instructive to compare the contributions of individual geometries in 
the two approaches. The physically important quantity is i?3(w), and one would not 
necessarily expect i?3a(w) and R^bi^j) to be the same in the two approaches. However 
it is of interest to see where the difference between the two approaches comes about. 
Firstly we will consider the impurity diagram result for i?3a(w), this result is the 
contribution of the path geometry shown in Fig. 2c, 
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In (2 + e)-dimensions where e is small, this becomes, 
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Comparing (45) and (54), we see that the impurity diagram approach has an extra 
term not present in the semiclassical result, just as in the previous order. Like in the 
two loop case the semiclassics gets the wrong prefactor on the term it does reproduce, 
in this case a factor of —3/4. However much worse than this, the second term in 
(54) has a finite real part in two dimensions, unlike the second term in the two loop 
diagrammatic result (41). This means the semiclassical result _R|^ is very different 
from the diagrammatic result for R^a- 

We now consider the impurity diagram result for Ry,i,{ijj) in dimension d = 2 + e, 
where e is small. We find, 
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where Ii{uj) was evaluated using (39), g is the dimensionless conductance given by 
(33), and s — uj/A. In the limit that e ^ it is clear that (56) is not zero, therefore 
we find that the semiclassical and impurity diagram technique results are completely 
different. 

The difference between the results of the two approaches can be summarised 
as follows. In the impurity diagram technique there is no (ln|s|/|s|) term because 
the terms in (55) and (56) which are of order (e~^(— is)^'^/^"-'^) exactly cancel each 
other. In the semiclassics there is no cancellation, and therefore the presence of a 
prefactor of order (e^^) means one must expand ((— zs)"^"^/^^^) to order (e^) to get 
the e-independent term. It is this that generates the (ln|s|/|s|) term. 
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8. The weak diagonal expansion for a quantum chaotic system in the 
ergodic regime. 

It is straightforward to apply the rcsuhs of the weak diagonal expansion to a system 
in the ergodic regime. If a system is ergodic, a particle is equally likely to be found 
anywhere in the system's phase space. This means that the propagation probability 
for a particle in such a system is featureless and hence the same for all ergodic systems. 
The method of evaluating the weak diagonal expansion presented here is cast purely in 
terms of this propagation probability and therefore predicts the same level statistics 
for all chaotic and disordered systems in the ergodic regime. Given this statement 
one would like to know two things: (i) are the results of the method discussed in this 
paper to be believed when applied to the ergodic regime, (ii) are the level statistics 
predicted by this method the same as those of RMT? A partial answer to the first 
question would be that it would be difficult to believe the results of the method in 
the ergodic regime, considering the failure of the method in the diffusive regime of 
a disordered system. However we can present a better answer to the first question 
by considering the answer to the second question. While the conjecture that chaotic 
systems have RMT level statistics [9] is generally considered unproven, this is not the 
case for disordered systems, for which a proof exists [8]. Therefore if this method's 
prediction for the ergodic regime of a disordered system does not agree with RMT 
[10], then this can be considered a proof that the method does not work in the ergodic 
regime. 

To evaluate the weak diagonal approximation results in the ergodic regime, 
we simply note that in this regime, the average propagation probability, 
{Je (f ', p', t\ f , p)) is just the inverse volume of the phase space surface with energy, E. 
Thus for a system of volume L , containing scmiclassical scatterers, the propagation 
probability is, 

(/^(r',p',i;r,p)) = — 1^ (57) 

We can now construct an expansion for the SFF, K{t), in the small parameter t/tn. 
This expansion can then be compared to the small t expansion of the SFF for random 
matrices [10]. If each term in the weak diagonal expansion does not agree with the 
corresponding term in the expansion for random matrices, then we have proven that 
this method of evaluating the weak diagonal approximation is wrong. 

The first term in the SFF is linear in t/tn and is given by the strong diagonal 
approximation which agrees with the leading order RMT result [14]. It is well known 
that for /3 7^ 2 RMT predicts corrections in powers of t/tn, while for /3 = 2 the 
linear in t behaviour is exact for t < tn. The weak diagonal expansion, if correct, 
must reproduce the higher-order contributions for [3—1 and yet yield no higher-order 
contributions for j3 = 2. 

In the ergodic regime with time-reversal symmetry {(3—1) the leading weak- 
diagonal contribution is the two loop term given by substituting (57) into (23). 
Comparing this result with the (3—1 RMT result, we immediately sec it is positive, 
while the RMT result is negative. Hence we find that this method of evaluating the 
weak diagonal approximation predicts level statistics different from those of GOE 
random matrices. 

In the (3 = 2 case, the leading order weak diagonal corrections are given by 
the three loop terms. It is well known from the impurity diagram technique that in 
disordered systems such contributions are individually non-zero, but that at each order 
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in the perturbation expansion they cancel each other. However when we calculate the 
three loop terms from the semiclassical approach, we find that both K^a and Ki,}^ (Fig. 
lc,d) are positive so there is no possibility of them cancelling each other. This means 
that again this method fails to reproduce the RMT results. 

We reiterate that, since we know that disordered systems have random matrix 
level statistics, the results of this evaluation of the weak diagonal approximation must 
be wrong. 

9. Concluding remarks. 

The main aim of this paper has been to cast the weak diagonal expansion in terms of 
semiclassics so that it could be applied to a generic chaotic system. If this could be 
achieved, the perturbation expansion would be a systematic way of finding corrections 
to the Berry diagonal approximation; thus in disordered systems such corrections 
give rise to weak localisation phenomena. If level statistics of chaotic systems are 
to be RMT, these corrections must be present to describe the difference between 
GOE statistics and the Berry diagonal approximation for times much shorter than 
the Heiscnberg time. 

Having presented this general scheme, we then developed a method of evaluating 
the terms in the weak diagonal expansion based on the Gutzwiller trace formula. 
Unfortunately, this method could not correctly reproduce the Hikami boxes which 
are regions in the phase space where quantum scattering is important. Thus it could 
neither reproduce the weak localisation expansion for disordered systems, nor give the 
corrections to the Berry diagonal approximation predicted by RMT. We do not believe 
that any existing semiclassical approach is capable of correctly evaluating the Hikami 
boxes for a generic chaotic system. A method exists [6] in which Hikami boxes were 
calculated for the system of semiclassical scatterers considered here. The problem 
with this method is that there is no obvious way to generalise it to generic chaotic 
systems. 

Finally we wish to encourage the use of disordered systems as a test bed for 
theories of quantum chaos. It is one of the few examples of chaotic systems whose 
quantum behaviour is well understood. In particular, any theory which claims to 
explain the RMT behaviour of quantum chaotic systems should be tested by checking 
if it produces the correct level statistics for disordered systems in the diffusive regime. 
It also appears that it is not enough for a theory to get the correct weak localisation 
correction to the conductivity of a diffusive system - this is too simple a test to be 
trusted. A complete semiclassical theory must be able to correctly reproduce higher 
order corrections to the conductivity and the TLCF. 
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